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We present first-principles simulations of single grain boundary reflectivity of electrons in noble 
metals, Cu and Ag. We examine twin and non-twin grain boundaries using non-equilibrium Green's 
function and first principles methods. We also investigate the determinants of reflectivity in grain 
boundaries by modeling atomic vacancies, disorder, and orientation and hnd that both the change 
in grain orientation and disorder in the boundary itself contribute significantly to refiectivity. We 
find that grain boundary reflectivity may vary widely depending on the grain boundary structure, 
consistent with published experimental results. Finally, we examine the refiectivity from multiple 
grain boundaries and find that grain boundary refiectivity may depend on neighboring grain bound- 
aries. This study raises some potential limitations in the independent grain boundary assumptions 
of the Mayadas-Shatzkes model. 
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I. INTRODUCTION 

As semiconductor technology continues scaling, metal 
interconnects must scale with minimum feature size (25 
nm in the generation currently in development) to con- 
nect to nanometer scale devices. However, in nanoscale 
metal wires, conductivity can degrade by more than 50% 
compared to the bulk metal [iTIll- As a result, both 
circuit delay times and power consumption may soon be 
dominated by interconnect [2]- 

The causes of this degradation are believed to be grain 
boundary, surface roughness, and interface scattering 
[J, As interconnect cross sections decrease, the inter- 
actions at the surface or interface with other materials 
become more important. Also, grain size is observed to 
scale roughly with wire thickness 0, H^, HJ]. However, 
there is limited understanding of the relative importance 
of these scattering sources on the observed degradation 
[l[ , nor is it known to what degree defects and impurities 
modulate these effects. Widely used semi-empirical mod- 
els of surface [1, 01 and grain boundary scattering |5i] and 
quantum models of surface scattering Q have been de- 
veloped, but there has not been sufficient understanding 
of grain boundary scattering at the atomic level. Else- 
where, two of us have studied the effects of wire surface 
roughness [l| and barrier/adhesion/seed layer scattering 
on conduction. 

In this paper, we present atomic-scale modeling and 
analysis of grain boundary reflectivity in metals and com- 
pare these results with experimental data. We also an- 
alyze the structure of grain boundaries, investigating in 
detail the determining factors of scattering. The organi- 
zation of the paper is as follows. First we present a sum- 



mary of our method and our results for single-boundary 
reflectivity. Next, we give comparisons to published re- 
flectivity measurements, both for multiple- and single- 
boundary measurements. Then we give detailed analy- 
ses of the causes of grain boundary reflectivity, breaking 
these into effects of a mismatch between Bloch states 
in the two crystallites and the effects of disorder and 
defects. Finally, we raise some concerns regarding the 
standard model in the literature for relating microscopi- 
cally calculated reflectivities with macroscopic resistivity, 
and present our investigations into this model's validity, 
as well as proposing some suggestions on how the model 
could be improved. 



II. METHOD 

We use the Non-Equilibrium Green's Function 
(NEGF) method and the Landauer formalism in this 
work p\. We use the code Atomistix[22'| to implement 
the calculations. In the Landauer formulation the re- 
sistance in perfect metal crystals at absolute zero (bal- 
listic conductors) is contact resistance (the reciprocal of 
Sharvin conductance). This is caused by the flnite num- 
ber of transverse modes per unit area with cutoff less 
than the Fermi energy. The number of modes is equiva- 
lent to the cross sectional area of the Fermi surface in a 
plane normal to the transmission direction: 
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where Gs is Sharvin conductance, Tg is ballistic trans- 
mission, ri±_ oc S/kE is a unit vector normal to the Fermi 
surface, z is the transmission direction, M{Ep) is the 
number of forward-moving modes with E = Ep , and 
the integration domain is the set of points on the Fermi 
surface with ri^ ■ z > 0. 



2 




III. COMPARISON TO EXPERIMENT 



In this section, we detail the agreement of our calcu- 
lated grain boundary reflectivity with experimental re- 
sults. The difference in length scales (along the transport 
direction) between those systems that can currently be 
simulated and actual wires necessarily introduces some 
uncertainty into the comparison. Here we compare re- 
sults despite this theory / measurement gap, but we re- 
turn to the issue in section IVl 



FIG. 1: Structure of (a) twin (210) / (120) and (b) twin 
(320) / (230) boundaries in Cu and (c) non-twin (100) / (110) 
boundary in Ag. Transmission is from left to right in all 
figures. 



TABLE I: Summary of simulation results. All systems are 
periodic in plane normal to transmission direction. 



System Relaxed? _Rcu (%) ^Ag ( 


Twin (210)/(120) 


Y 


17 


12 


Twin (320)/(230) 


Y 


13 


14 


Non-Twin (111)/(110) 


Y 


47 


36 


Non-Twin (110)/ (100) 


Y 




46 


Non-Twin (111)/(100) 


Y 


19 


16 


Vacancy (39.2 A^)"^ 


N 


8 




Vacancy (19.6 A^)"^ 


N 


16 





In this paper we investigate the reflection probabil- 
ity R for various structures. In our notation, the pres- 
ence of scattering modifies transmission from its ballistic 
value Ts ( Eq. (P) ) to T = Gh/2e^ = (I - R) M{Ef). 
We use computer simulations based on density functional 
theory^ within the local density approximation (LDA) 
[lol | and NEGF to estimate electron transmission T at 
K. For structural relaxation, we use the total energy pseu- 
dopotential method and perform the relaxation within 
LDA, using augmented wa ve p seudopotentials with pe- 
riodic boundary conditions! Ill [T2j. Due to the periodic 
boundary conditions, we use at least 10 atomic layers 
for both sides of the grain to avoid the effect of repeat- 
ing images on the configurations of grain boundaries. For 
controlled atomic defect studies due to vacancies and dis- 
order, we do not relax the structures. 

We simulate twin (coincidence site lattice, or CSL) and 
non-twin grain boundaries in two FCC (face centered cu- 
bic) metals, Cu and Ag. We prepare twin boundaries 
for the two angles with smallest supercells, correspond- 
ing to (210) / (120) and (320) / (230), as shown in figures 
[IJa) and (b). We also prepare non-twin boundaries for 
(111) / (110), (110) / (100), and (111) / (100) as demon- 
strated in figure [ijc). The refiectivity simulation results 
are summarized in table ID 



A. Wire Resistivity Measurements 

In the literature, the Mayadas-Shatzkes (MS) model Q 
is the most widely used analytical model to extract grain 
boundary reflectivity from experimental measurements 
of resistivity p. The MS model describes metal conduc- 
tivity as a function of boundary reflectivity R and grain 
size D. The model's basic assumptions are that grain 
boundaries are randomly distributed, that all reflections 
are specular and occur with probability R independent 
of incoming momentum, and that transport is semiclas- 
sical between boundaries. As we will describe in section 
rVl several considerations lead us to question the valid- 
ity of these assumptions. Still, we compare our results 
for R to measurements of resistivity p by using MS, as 
it is currently the standard model of p as a function of 
microscopic properties. 

The MS model agrees with an even simpler model 
we constructed based only on the assumptions of the 
Landauer formula [8] and multiple scattering from grain 
boundaries with average reflection probability R. This 
is important because it allows us to extend the model's 
predictions to more general cases. For resistivity of a 
conductor with both bulk scattering of mean free path 
At, and grain boundary scattering, both theories give the 
resistivity augmentation over bulk as: 
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in our simplifled model and k « 1.39 in MS [l3|. Thus, 
at room temperature and _D = 45 nm, a grain boundary 
refiectivity of 20% increases resistivity by ^31% over the 
bulk value. 

We use the MS model to compare our results for R at 
K to experimental results at 5 K, in the regime where 
bulk scattering effects are minimal. Although Cu is more 
important for integrated circuits, more experimental data 



is available on Ag. The low temperature experiments [13 1 
indicate R w 25% for Ag, in comparison to our values of 
R ~ 12% for twin boundaries and R from 16% to 46% 
for non-twin boundaries (table |l|. 
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A survey of experimental results indicates R values 
for Cu and Ag (computed with MS) in the range of 24% 
to 46%, in good agreement with our predictions. 



B. Single-Boundary Resistance Measurements 

A more direct comparison can be provided by a few ex- 
periments measuring grain boundary resistance directly 
in a metal. Schneider et al.^VJ\ measured potential dif- 
ference across single grains in Au (an fee noble metal like 
Cu and Ag), and found reflectivity in the range 70% to 
90%, depending on orientation. NakamichifTsI measured 
interface resistivity pgb = RgbA (where Rgh is interface 
resistance) in Al for various single grain boundaries. To 
compute reflectivity from Nakamichi's results, we esti- 
mate the ballistic conductance per unit area Gc/A for Al 
by Eq. ([T]) assuming a spherical Fermi surface. We then 
compute reflectivity with 



A/Gc 



{Pgb + A/Gc) 



(3) 



Aggregating Nakamichi's results and analyzing according 
to this expression, we find most of the measured twin 
boundary results in the range R = 0% to 27% and non- 
twin results in the range R — 36% to 51%. This range 
of reflectivity variations agrees well with our findings, 
particularly given the different materials used. 



IV. DETERMINANTS OF REFLECTIVITY 

Although real grain boundaries in metals are not neces- 
sarily two-dimensional plane defects nor form orthogonal 
interfaces to the transmission direction, we have isolated 
the effects and simulated them individually. 

Since a grain boundary is the interface between two 
crystallites, there are two broad categories of scattering 
that could occur: 1) scattering caused by the misaligned 
crystal orientation of the two grains, a category we re- 
fer to as orientation effects] and 2) the atomic structure 
in the interface itself, which we call atomic position ef- 
fects. We present simulations and analytical arguments 
to quantify these two effects. 



A. Orientation effects 

We wish to understand the contribution to scattering 
from the change in orientation across grains, Rs. An ideal 
coincidence site lattice (CSL) twin boundary has zero 
thickness, so the only possible cause of scattering is the 
abrupt change in orientation. We therefore consider our 
results for reflectivity of unrelaxed twin CSLs, R ^ 15%, 
an estimate of the orientation effect. 

In table [H most of the non-twin reflectivities are higher 
than the twin ones. This is likely due partly to atomic 



position effects (since the interface is less sharp). But 
it may also be due to the difference in Bloch bases in 
the two grains, particularly if incident states are poorly 
approximated by states in the other grain. 

We therefore simulate transmission in several differ- 
ent grain orientations. We estimate the orientation ef- 
fect, due to impedance mismatch between different sets 
of Bloch states in the two crystallites, by the relative 
difference 

\AT\/T 

in transmission across the boundary. (This is by analogy 
to a simple ID quantum potential step. 



R 



k-k' 
k + k' 



where k and k' are the momenta of the states on either 
side of the step.) We find that transmission in the (100) 
and (111) directions are similar, while transmission in the 
(110) direction is significantly higher. We find transmis- 
sion for two intermediate orientations (210) and (320) 
to be between that of (100) and (110). This shows a 
large dependence of T on orientation for a perfect crys- 
tal, although in a realistic system, we expect that bulk 
scattering would diminish the orientation effect. 

To confirm this large dependence of transmission T on 
orientation we compare the simulation results with an 
analytical estimate of T as a function of orientation in 
Cu. The integral in Eq. ([IJ may be evaluated numer- 
ically for different directions z. The deviations from a 
spherical Fermi surface then give the effect of grain ori- 
entation on ballistic transmission. In the Fermi surfaces 
of Cu and Ag, the eight (111) directions each contain 
a "neck" that intersects the Brillouin Zone edge. These 
necks contribute to Eq. (P) for (110) transmission. Such a 
calculation has been carried out by Xu et al. [l6j , giving 
T(110)/T(100) = 1.07 and T(lll) « T(IOO). This quali- 
tatively agrees with our NEGF-based prediction that the 
transmission in (111) is similar to (100) and greater in 
(110) relative to (100). 



B. Atomic Position Effects 

We isolate the atomic position effect in the interface 
by simulating structures with specific structures and ex- 
tracting the magnitude of the reflectivity, Ra. Our initial 
assumption is that total reflectivity from atomic position 
in the boundary can be decomposed as a sum of reflec- 
tivities from individual defects (individual moved atoms, 
etc.). 



Ra — ^ ^ Ra.ij 



(4) 



when the individual reflectivities Ra,i are sufficiently 
small. 
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FIG. 2: System with plane of vacancies at density (39 A^)~^. 
Note that structures are periodic out of the page, but alter- 
nate atoms are inequivalent. 




O Simulation, RMS displacement 0.7 A 
n /(n + 8 81| 

* Simulation, RMS displacement 0.24 A 
n /(n +5D.0I 




Number n of disordered iayers 



FIG. 3: System with 6 disordered layers (layers 3-8 of those 
shown) . 



FIG. 4: Total reflection probability R vs. number n of dis- 
ordered layers. Systems with smaller n have a subset of the 
disordered region present in systems with larger n. There is 
significant variability evident due to magnitude of displace- 
ment. 



1. Vacancies 

We model vacancies in grain boundaries by simulating 
structures with a single interface containing vacancies as 
shown in figure [5J We examine two different vacancy 
densities, one or two vacancies per 39 (3 unit cells) of 
cross-sectional area. Results are R = 8% and 16% (table 
HI, scaling linearly with vacancy density consistent with 
additive reflectivities. The scattering cross section is on 
the order of magnitude of the area of the missing atoms, 
as one missing atom per 39 corresponds roughly to 
one vacancy in 12 atoms (3 unit cells) of the interface 
giving an analytical expectation of i? ~ 1/12 = 8.3%, in 
agreement with the calculations. 



magnitude of the atomic position effect: 
Ra - (0.08) (39 A'), 

with del the defect area density in the boundary, and 
5% < Ra < 30%, depending on the magnitude of disorder 
in the boundary. The objective of these estimates is to 
approximately predict Ra, as function of structure. Both 
estimates reflect a strong dependence of grain boundary 
reflectivity on atomic position (e.g. gaps, relaxation). 



V. VALIDITY OF MAYADAS-SHATZKES 
MODEL 



2. Disorder 

To further investigate the effect of atomic position on 
reflectivity, we simulate layers of disordered Cu atoms as 
shown in figure [3] to isolate the impact of crystalline or- 
der on conduction. The disordered atoms are displaced 
by normally distributed random vectors with root mean 
square (RMS) mag nitude 0.24 A or 0.70 A. We change 
the number n of such layers and expect from theory 
that R = n/{n + no) with a constant ng, since localiza- 
tion length is long compared to our system sizeQ. We 
summarize the results in figure |4l These disordered re- 
gion simulations give insight into the effect of non-linear 
boundaries on refiectivity and show that the impact of 
non-lattice site atomic positions is significant. 

These structures are not the only possible departures 
from crystalline order, but provide an estimate of the 



We also attempt to check the validity of the Mayadas- 
Shatzkes (MS) model assumptions using model systems 
with multiple grain boundaries. The aim of this inves- 
tigation is to raise the question: are the assumptions in 
the MS model, such as the characterization of all grain 
boundary scattering in a metal by a single reflectivity pa- 
rameter, tenable? This undertaking is motivated partly 
by a lack of consistency in published reflectivities for the 
same metals. 

The Mayadas-Shatzkes model makes many assump- 
tions, which we summarize here: 

• Transport is semiclassical between grain bound- 
aries. 

• Grain boundaries are perpendicular to transport. 

• Grain boundaries are translation-invariant along 
the boundary. 
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• Transmission at grain boundaries can be character- 
ized by a single parameter R. 

• All boundaries in a sample are identical (same R). 

Some of these assumptions are clearly objectionable, and 
it seems likely that some may affect the final result. 
For example, consider the assumption of translation- 
invariant boundaries. In the presence of bulk scattering, 
velocities are randomized within a few mean free paths 
of the boundary, so the assumption of specular reflection 
or undeflected transmission may not matter, as long as R 
matches the average reflection probability. But if grain 
size is comparable to mean free path (hard to avoid with 
wire thickness ^ 25 nm < room temperature mean free 
path of 39 nm), or several boundaries are particularly 
close together, then the particular states transmitting 
through one boundary and incident on a second may af- 
fect overall transmission. We simulate just such a system 
here. 

We did simulations with two (210) / (120) twin bound- 
aries separated by from 2.4 nm to 4.6 nm (figure [5]). Al- 
though this differs markedly from an infinite system as 
considered in the MS model, the discussion surrounding 
Eq. ^ justifies comparing the results using our simplified 
Landauer transport model for multiple reflection between 
just two identical boundaries. This model gives the total 
reflectivity as 



Rt ^2R/{1 + R), 



(5) 



where R is the reflectivity of a single boundary[8] (assum- 
ing averaging over angles of reflection destroys the co- 
herence in the electron wavefunction between the bound- 
aries). Using our single-boundary results (table H]), we 
anticipate Rt = 29% for Cu and Rt = 21% for Ag. 

We find our simulation results differ from the smallest 
to largest grain size as shown in table |lll We attribute 
this change in reflectivity to interactions between bound- 
aries. We also find the simulated Rt is somewhat lower 
than that predicted by (O . This may be explained by the 
observation that transmission depends strongly on initial 
momentum in our single-boundary simulations. The first 
boundary may act as a filter, letting only those states 
with highest transmission through to the second bound- 
ary. If a real sample contained a pair of similar bound- 
aries with spacing small compared to the bulk mean free 
path (like our model systems) , one might therefore expect 
the MS predictions to fail. The MS model also fails to 
consider nonspecular reflection and transmission of elec- 
trons. 

Still, this explanation would be affected by bulk scat- 
tering or variety in boundary type, factors which render 
the semi-classical MS treatment more tenable. The fail- 
ure of MS in the system here could be criticized on the 
grounds that the simulated system is unrealistic. Our re- 
sponse is that small grain separation is becoming likelier 
as wire dimensions shrink, while our neglect of phonon 
scattering is rendered more realistic by small grain size 
and the inclusion of disorder in our simulated boundary. 



FIG. 5: System used for validation of MS model. Two twin 
grain boundaries are separated by an ideal (ballistic) crys- 
talline region. 



TABLE II: Total grain reflectivity of two twin grain bound- 
aries with various grain sizes (in units of lattice constant 
acu ~ 3.6 A, SAg = 4.1 A) for Cu and Ag. 

Grain size/a Rt for Cu (%) Rt for Ag (%) 



6.71 
8.94 
11.16 



23.6 
25.8 
24.5 



19.3 
14.7 
17.4 



Expected 



29 



21 



The ideal test of MS would be a first-principles simula- 
tion with multiple, different grain boundaries separated 
by a more realistic grain size, together with bulk scat- 
tering, but this unfortunately is not yet practical. The 
goal of our test has simply been to narrow the theory- 
measurement gap with the computing resources currently 
available. Based on these analytical and computational 
considerations, we propose that an MS-like model with 1) 
multiple reflectivity parameters for different incident mo- 
menta, 2) a treatment of the statistical likelihood of bal- 
listic propagation between boundaries, and 3) allowance 
for deflected transport and grain boundary angles would 
capture more relevant physics and probably give more 
consistency across measurements. 



VI. DISCUSSION 

Our findings show that grain boundary scattering is a 
significant source of resistivity for sufficiently small grain 
size. However, the damascene process can deposit inter- 
connects with average grain size of order the minimum 
feature size or larger |18l42l| . so grain boundary scatter- 
ing can be reduced but not eliminated. Grain boundary 
scattering is still probably more important than the rela- 
tively small contribution from surface roughness scatter- 
ing [1] , and even the moderate contribution from barrier 
layer interface scattering \Jj . 

The damascene process is capable of depositing Cu in- 
terconnects with grain size larger than the line thickness. 
For example, Geiss and Read 18j report an average grain 
diameter of 315 nm for 100 nm damascene Cu lines. Sim- 
ilarly, Paik et al.\l^ measured grain sizes in the range 
125 - 275 nm for ~170 nm line thickness. Both Carreau 
et aZ. [20*1 and Steinhogl et al. [2l| measured somewhat 
smaller grains for thinner damascene-deposited intercon- 
nects, and observed that grain size does indeed scale with 
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thickness for the thinnest wires. 

These results suggest that grain sizes on the order 
of the minimum feature size are readily achievable, and 
that annealing and overgrowth have the potential to give 
larger grain sizes. If we assume that grain size is equal 
to minimum feature size, Eq. ([2]) gives a resistivity aug- 
mentation of 42% for a minimum feature size of 32 nm 
and average R — 20%, and augmentation of 110% if the 
average of R is 40%. On the other hand, if a grain size 
of 3 times the minimum feature size can be achieved (as 
suggested by the results of Geiss and Read) then resis- 
tivity augmentation is only 14% for R = 20%, and 38% 
for R = 40%. 



VII. SUMMARY 

We have presented the first simulations of reflectiv- 
ity for relaxed twin and non-twin grain boundaries in 
Cu and Ag. Our results agree with the experimental re- 
ported range of reflectivity [isl . [l3 and with individual 
grain boundary measurements |l5l . [l7j . To gain insight 
into the mechanisms of grain boundary reflectivity and 
the impact of the non-planarity of real grain boundaries. 



we also investigated the effect of vacancies, orientation, 
and disorder. We observed that all three contribute sig- 
nificantly to reflectivity. Our predicted dependence of 
reflectivity on grain boundary type and isolated vacancy, 
orientation, and disorder effects may explain the wide 
range of variations in the experimentally data. In prob- 
ing the utility and extendibility of the Mayadas-Shatzkes 
model we found that the assumption of a one-parameter 
reflectivity averaged over all grain boundaries and ini- 
tial states failed to accurately estimate reflectivity from 
closely-spaced multiple grain boundaries. Improvements 
in the analytical models to account for deviations from 
additivity of multiple boundaries and the impact of grain 
boundary type and non-planarity are needed. Based on 
our studies, it is clear that larger-scale rigorous quan- 
tum models are needed to capture more realistic line and 
grain boundary structures. 
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